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The dynamics of a two-qubit system is considered with the aim of a general categorization of the 
different ways in which entanglement can disappear in the course of the evolution, e.g., entanglement 
sudden death. The dynamics is described by the function n(t), where n is the 15-dimensional 
polarization vector. This representation is particularly useful because the components of n are 
direct physical observables, there is a meaningful notion of orthogonality, and the concurrence C 
can be computed for any point in the space. We analyze the topology of the space S of separable 
states (those having C = 0), and the often lower-dimensional linear dynamical subspace D that 
is characteristic of a specific physical model. This allows us to give a rigorous characterization 
of the four possible kinds of entanglement evolution. Which evolution is realized depends on the 
dimensionality of D and of D n S, the position of the asymptotic point of the evolution, and 
whether or not the evolution is "distance-Markovian" , a notion we define. We give several examples 
to illustrate the general principles, and to give a method to compute critical points. We construct 
a model that shows all four behaviors. 
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I. INTRODUCTION 



Entanglement is one of the most intriguing aspects of 
quantum physics and is known to be a useful resource for 
quantum computation and communication [l|, |2[ . How- 
ever, its structure and evolution in time are not fully 
understood even for simple systems such as two qubits 
where a relatively computable entanglement measure, the 
Wootter's concurrence C, is available Q. The difficulty 
resides in the high dimensionality of state spaces (fifteen 
for two qubits) and the non-analyticity of the definition 
of the entanglement measure C. 

In the presence of external noise, pure states become 
mixed and entanglement degrades. These are distinct 
but related issues. The completely mixed state (density 
matrix proportional to the unit matrix) is separable: it 
has C = 0. Pure states, on the other hand, can have 
any value of C between and 1 inclusive. Purity (mea- 
sured, for example, by the von Neumann entropy) tends 
to decrease monotonically and smoothly with time under 
Markovian evolution. The same is not true for entan- 
glement. Apart from the expected smooth half-life (HL) 
decaying behavior, the sudden disappearance of entan- 
glement has been theoretically predicted and experimen- 
tally observed 0, [3 • Widely known as entanglement sud- 
den death (ESD), this non-analytic behavior has been 
shown to be a generic feature of multipartite quantum 
systems regardless of whether the environment is mod- 
eled as quantum or classical [a-[l9|. While the mono- 
tonic decrease of C(t) is usually associated with Marko- 
vian evolution, non-Markovian evolution can also lead to 
entanglement sudden birth (ESB). It is believed to be 
related to the memory effect of the (non-Markovian) en- 
vironment [fj|, flSl — l2l| . Although most investigations have 
been focused on two-qubit systems and we will also focus 
on this case, ESD and ESB have been shown to exist in 
multi-qubit systems, and even in quantum systems with 



infinite dimensional Hilbert spaces, such as harmonic os- 
cillators fUS 



Our aim in this paper is to formulate the problem of 
the evolution of entanglement in two-qubit systems in the 
polarization vector space, and to show that this formula- 
tion leads naturally to a categorization of entanglement 
evolutions into four distinct types, generalizing and mak- 
ing precise the concepts of HL, ESD, and ESB behaviors. 
It turns out that these categories arc consequences of cer- 
tain topological characteristics of a model. To show this, 
we proceed as follows. Sec. |H] characterizes in detail two 
manifolds in the polarization vector space: the manifold 
of admissible physical states and the manifold of separa- 
ble states. Sec. IIIII presents the concept of a dynamical 
subspace - a manifold that is associated with a physical 
model, and then gives several concrete examples of mod- 
els of increasing complexity. In Sec. IHII we also show how 
to compute critical points: parameter values that sepa- 
rate one behavior from another. In Sec. IIV[ we prove 
that our categorization scheme is exhaustive. Sec. |V| 
presents the final results and discussion. 



II. ENTANGLEMENT IN POLARIZATION 
VECTOR SPACE 



The universality of the various entanglement behav- 
iors suggests that they are derived from some structural 
property of entanglement in the physical state space, and 
that the system dynamics can be viewed as a probe of 
that property. To state this property precisely, we need 
to first characterize the space of all admissible density 
matrices p, or equivalently, the space of all admissible 
polarization vectors. 

For two qubits, the polarization vector n = 



(m, ri2, . . . , TI15) is defined by the equation 



P=- A h 



1 15 



(1) 



where I4 is the 4x4 unit matrix and the pi are the 
generators of SU(4), satisfying 



M« 



4, 



Tr Hi = 0, Tr /i,^- = 4(5; 



(2) 



For our purposes the pi are most conveniently chosen as 

Pa/3 = <y a ®<7fi, (3) 

where a a acts on the first qubit and er^ acts on the sec- 
ond qubit. a and /? sum over the 2x2 unit matrix / 
and the Pauli matrices X, Y and Z. Thus, in Eqs. Q] 
and [21 i is regarded as a composite index of a and /3, but 
the a a = Ofi = I term is singled out. This space has the 
usual Euclidean inner product (which corresponds to the 
Hilbert-Schmidt inner product on the the density matri- 
ces) , and the inner product induces a metric and a topol- 
ogy in the usual fashion. The components of ft are phys- 
ical obscrvables and can be calculated by rii = Tr ppi. 
For example, the average value of the z-componcnt of the 
spin of the first qubit is (Z <X> I) = Tr [p(Z <X> /)] = nzi- 
The six components nix, nrr> niz, rixi, nyi and nzi 
represent physical polarizations of spin qubits. The other 
nine components (rixx, etc.) are inter-qubit correlation 
functions. The most common name for ft is "polarization 
vector", but "coherence vector" and "generalized Bloch 
vector" are also in use. We note that different normaliza- 
tions in Eqs. Q] and [2] are used in the literature |29l - l33| . 
The generators jXi satisfy 



HiHj 



= SijI- 



(ifijk + di jk )pk 



(4) 



where fijk is totally anti-symmetric and dijk is totally 
symmetric. These structure constants can be found in 
Appendix [A] 

Eq. [1] holds for a 4-level system. It has an obvious 
generalization to iV-level systems; the pi just become 
the N 2 - 1 generators of SU(N). For N = 2, ft is the 
usual Bloch vector in a real 3-dimcnsional vector space. 
It is important to stress that the correspondence between 
p and ft is one-to-one; they give completely equivalent 
descriptions of the physical system. Certain physical 
concepts have geometric interpretations when stated in 
terms of ft; as we shall see below. This is not so true of 
p. In our opinion, ft is the more convenient quantity for 
most purposes, p has been the traditional language in 
which to describe mixed states, but some experimental 
groups now favor ft [3J, [35| ■ 

We shall refer to the set of all admissible ft as M, the 
state space. What shape does M have? Eq. Q]guarantees 
that p is Hcrmitian and has unit trace. To guarantee that 
p is positive (all its eigenvalues are non- negative) , we also 
need the condition that all coefficients a,- of the charac- 



/v 



tcristic polynomial det(x In — p) = J2j=o(~^) 
are non- negative [36|. Note ciq = 1 by definition. 



djX 



N-j 



For two-qubit systems there are four of them, which 
are 



l!oi = Tr p = 1, 

2!a 2 = 1 - Tr p 2 , 

3!a 3 = 1 - 3Tr p 2 + 2Tr p 3 



(5) 
(6) 

(7) 
4!a 4 = 1 - 6Tr p 2 + 8Tr p 3 + 3 (Tr p 2 ) 2 - 6Tr p 4 . (8) 

Note that a\ > is trivially satisfied for all density ma- 
trices. 

For one qubit N = 2, ft is the usual Bloch vector, and 
only the a 2 > constraint applies. Thus the positivity 
requirement is that \n\ < 1 and M is the familiar 3- 
dimensional spherical volume. For the 2-qubit N = 4 
case that we are concerned with, there arc cubic and 
quadratic inequalities to be satisfied, so the surface that 
bounds M is not so simple. The main point, however, is 
that M is convex: the line joining any two points in M 
is also in M. This follows from the convexity argument 
for p: if p\ and pi are positive, then so is sp\ + (1 — s)p% 
for all < s < 1. This argument clearly also holds for ft. 

All of the positivity requirements can be written in 
terms of ft, but the higher-order ones arc fairly compli- 
cated. The requirement 1— Tr p 2 > is of particular 
interest, since it has a simple expression in terms of ft : 



< 1 - Tr 



1 

16 



P 



15 15 

Tr I4, + 2 y^ njTiHi + 2J n.injTrpiPj 



= 4 - 4 I"! 2 ' or I"! 2 - 3 - 

Hence the vectors in M lie within a sphere of radius s/3. 
Technically, M is a 15-dimcnsional manifold with bound- 
ary. We will follow physics usage and also employ the 
term "space" for M , though of course it is not closed 
under vector addition. Note that pure states satisfy Tr 
p 2 = 1, so the pure states are a subset of the 14-sphere 
in M with \ft\ = v3- To be more specific, the two-qubit 
pure states |'0)(V'| ar c of measure zero on that sphere 
since they can be parametrized by 6 real parameters 

|V>) =cos6»i|00)+e^ 1 sin6'isine' 2 |01) 
+ e^ 2 sin6»icos6i2cos6l 3 |10) 



o «03 



sin B\ cos 02 sin 6*3 1 11) 



An overall phase has been dropped in writing \tp) since 
it does not appear in p. 

Further insight into the shape of M can be gained 
by noting that M must be invariant under local unitary 
transformations (rotations of one spin at a time) , which 
means that M has cylindrical symmetry around the 6 
singlc-qubit axes. This is verified by making some 2- 
dimensional sections of M with exactly two components 
of ft non-zero. In contrast to Ref. [37| where a different 
basis was used (generalized Gell-Mann matrices), we find 
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FIG. 1. (Color online) Two dimensional sections of the full 
two-qubit state space M, keeping only two components of the 
polarization vector non-zero. Only two inequivalent sections 
are possible: the disc section (D) and square section (S). The 
physical states are colored grey. a2 = is outlined as black, 
<23 = blue and en — red. All physical states in these two 
dimensional section are separable states. 



only two types of shapes, as shown in Fig. Q] and tabu- 
lated in Table [I] Using the structure constants dijk and 
fijk, it can be shown that the square and disc sections 



are the only possibilities along the n. 



plane. When 



Hi commutes with pj the section is a square and pi anti- 
commutes with pj the section is a circular disc. Details 
can be found in Appendix [X] 

The discs correspond to the local rotations between 
single-qubit-type axes, such as the nix ~ n iY section. 
If, on the other hand we rotate from a definite polariza- 
tion state of qubit 1 to a definite polarization state of 
qubit 2, we find a square cross-section; examples are the 
nix—nxi or nix — nyi sections. Rotations of ft that mix 
single-qubit-type and correlation-type directions can be 
of either shape; the nix — nxx section is square, while 
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TABLE I. Two dimensional sections of the full two-qubit state 
space M, keeping only two components of the polarization 
vector non-zero. Type S is the square section, while type D is 
the disc section. The table is symmetric thus the upper part 
is omitted. The diagonal entries do not correspond to two 
dimensional sections. 



the nix — nyy section is a disc. Finally, rotations be- 
tween correlation-type directions can have either shape. 
The nxx — ijy section corresponds to a local rotation 
of qubit 2; hence it is a disc. Rotations involving both 
qubits, such as that which generates the nxx — fiyy sec- 
tion, generally give square sections. 

We may conclude that M is a highly dimpled ball, 
perhaps most similar in shape to a golf ball. Its minimum 
radius is |n| = 1 and its maximum radius is |n| = V3- 

Since our aim is to quantify entanglement in M, we 
need an entanglement measure. We will employ C , the 
concurrence of Wootters Q . The concurrence varies from 
for separable states to 1 for maximally entangled state, 
i.e., the Bell-like states. It is defined as C = max{0, q}, 
and 

q = Ai — A 2 — A 3 — A4 , 

where A^ are the square roots of the eigenvalues of the 
matrix pabPab arranged in decreasing order and 



PAB = ((Ty®Vy)PAB(<7y 



(9) 



is a spin-flipped density matrix. p* AB is the complex con- 
jugate of the density matrix pab- It is not possible to 
write the function C (ft) in a simple explicit form unless 
further restrictions on ft apply [38l |. but it is clear from 
the form of the continuous function q (ft) and the pres- 
ence of the max function that C is a continuous but not 
an analytic function of ft. 

We next consider S, the manifold of separable states, 
which we define as those ft for which the concurrence 
vanishes: C(ft) =0. S is a subset of M and M\S is 
the set of entangled states. S includes the origin since 
q(6) = -1/2 and C(3) = max{0,q(0)} = 0. Since q(n) 
is continuous, S actually includes a ball of finite radius 



about the origin: it can be shown that if \n\ < l/yS, 
then ft £ S [23]. Thus the manifold of separable states 
has finite volume in M: S is also 15-dimcnsional. We 
will also refer below to the interiors and boundaries of 
M and 5 and denote these by Int(M), Bm, Int(S), and 
Bs- Since the various sets we encounter in this paper 
are not linear subspaces, we need the general topological 
definitions of "boundary" , "interior" and "dimension" . 
These may be found, for example, in Rcf. |40f . 

S is also a convex set. What else can we say about 
the shape of S7 It is easily seen that the surface of S, 
like that of M, is rather non-spherical. Indeed C = 
along any of the basis vector directions. Coupled with 
the fact that S is convex, we see that S must contain a 
large hyperpolygon with vertices at ft = (±1,0, 0, ...,0), 
n = (0, ±1, 0, 0, ..., 0), etc. C (ft) is invariant under local 
rotations, so it has the same hyper-cylindrical symmetry 
as M. Again we may consider 2-dimensional sections 
in order to understand the shape of the surface. Two 
examples are shown in Fig. [T] 

A simple-sampling Monte Carlo study shows that there 
are more entangled states than separable states in M. 
The details can be found in Appendix [B] 



III. DYNAMICAL EVOLUTION IN S 
A. introduction 

The dynamical evolution or trajectory of a quantum 
system is a function ft(t) with t G [0, 00) and n G M. 
The initial point is n (0) and, in the cases of interest 
here, the trajectory approaches a limiting point as t — > 00 
and we can define fioo = \w\t^oo n(t). The entanglement 
evolution is the associated function C (t) = C(n(t)). 
C G [0,1]. For studies of decoherence the main inter- 
est is in entanglement evolutions such that C(0) > and 
C(oo) = 0, i.e., the system starts in an entangled state 
and ends in a separable state. Four distinct categories 
of entanglement evolution of this type have been seen in 
model studies [42|. They are shown in Fig. O These 
four categories are topologically distinct, as may be seen 
by considering the set T = {t\C(t) = 0}. In category 
A, To is the null set; in category £, To is a single infinite 
interval; in category B, Tq is a set of discrete points; in 
category 0, To is a union of disjoint intervals. 

These categories also reflect how the trajectory tra- 
verses S and M\S. Entanglement evolutions in category 
A approach the boundary of separable and entangled re- 
gions asymptotically from the entangled side. The trajec- 
tories never hit S while the decrease in ent ang lement may 
or may not be monotonic, as seen in Ref. [431 ] . Entangle- 
ment evolutions in category B bounce off the surface of S 
at finite times but never enter S. Overall, entanglement 
diminishes nonmonotonically. Entanglement evolutions 
in category £ enter S at finite time and entanglement 
stays zero afterwards. This is the typical ESD behavior. 
Entanglement evolutions in category O give ESB: after 
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FIG. 2. Four categories of entanglement evolution. A: ap- 
proaching. B: bouncing. £ : entering. O: oscillating. 



ESD, entanglement suddenly appears after some dark pe- 
riod. 

We shall focus on models with associated linear maps 
A, i.e., p(t) = A(i)[p(0)]. More general non-linear models 
may be contemplated, but they seem to have unphysical 
features J44J. It is known that A(i) is completely positive 
(CP) if and only if there exists a set of operators {E a } 
such that [2[ 



A[p(0)]=Y,E a p(0)El 



(10) 



We require A to be trace preserving so that it maps den- 
sity matrix to density matrix. This condition is equiv- 
alent to the completeness condition ^ a E a E^ = I. In 
terms of the polarization vector, the dynamics is de- 
scribed by an affine map T(i) acting on the initial po- 
larization vector n(0), i.e. 



n(t) ee T(t) [n(0)] 

= T{t) n(0)+m(t) 



(11) 



where T(t) is a real matrix and fh(t) is a real vector 
[33J. rn{t) is zero for all time only when T(i) is unital, 
i.e., it maps to (in terms of A, the unital property 
means that A maps identity matrix to identity matrix). 
T(0) = I and m(0) = 0. 

Coherent dynamics is described by unitary transfor- 
mations on the density matrix (single Kraus operator). 
The dynamical map T(i) is then linear which translates 
to orthogonal transformations T acting on the polariza- 
tion vector n. Decoherent dynamics (multiple Kraus 
operators) is characterized by the nonorthogonality of 
the transfer matrix T . Markovian dynamics is conven- 
tionally defined by T possessing the semigroup property 
T(ii + £2) = T(t2)T(ii), which translates to 

T(t 1 + t 2 )=T(t 2 )T(t 1 ) (12) 

m(h + t 2 ) =T(t 2 )m(t 1 ) + m(t 2 ). (13) 

We shall adopt a slightly different definition of Marko- 
vianity for the present paper. An evolution will be said 



to be "distance Markovian" if \n(t) — n^] is a monoton- 
ically decreasing function. Note "distance Markovian" is 
a weaker condition than Markovian, though the two are 
usually equivalent. Given the semigroup property Eq. \T^ 
and Eq. [T3l we have 

\n{t) — fioo] = |[T(t)fJ(0) + m{t)} - [T(i)noo + rn(t)]\ 
^llTWHIn^-nool 
< |»f(0) — Tlool 

since all eigenvalues of T[t) have their norms in the range 
to 1, i.e., T(t) cannot increase the purity of the quan- 
tum state. 

Any model of an open quantum system defines a set 
of possible dynamical evolutions. This is done by speci- 
fying the equations of motion, which give T, fa and the 
initial conditions, which give n(0). We define the dynam- 
ical subspace D of a model as the set of all trajectories 
allowed by the set of initial conditions and the equations 
of motions. Eq. QT] shows that, as long as the set of all 
initial conditions is a linear space (the usual case), then 
D is a linear space intersected with M: we first choose a 
basis that spans the set of all possible n (0), then evolve 
this basis according to Eq. [TTJ giving a linear subspace 
in the space of all n. A precise and general definition of 
D is given in Appendix [Cj The set of admissible ft is 
then given by intersecting this linear subspace with M. 
D is a manifold of any dimension from 1 to 15 in the two- 
qubit case. We note dimD could be smaller than dimilf . 
This happens when both p(0) and E a (t) are expandable 
by identity and a true subalgebra a of su(N). dim!? 
is then equal to the number of independent elements in 
a. For example, if p(0) is a two-qubit "X-state" and the 
dynamics can be described by the action of Kraus oper- 
ators in the X-form, the dynamical subspace D will be 
7-dimensional [45J |. 

It is the nature of the intersection of D with S and 
the position of n^ relative to S that determines the cat- 
egories of entanglement evolution of a model. 

The aim of the remainder of this paper is to show how 
to determine the topological structure of D n S and the 
position Hoc for various illustrative models of increasing 
complexity, and then to deduce the possible entangle- 
ment evolutions from this information. We note that in 
general D can be determined without fully solving the 
dynamics. Thus it is possible to gain qualitative infor- 
mation of the entanglement evolution of the model with 
simple checks. 



B. Model D z 

Our first model consists of two qubits (A and B) with 
a Hciscnbcrg interaction and classical dephasing noise on 
one of the qubits. The Hamiltonian is 




FIG. 3. (Color online) Effective Bloch sphere representation 
for the dynamical subspace Da . The north and south pole are 
either {| 10) , |01}} or { 1 00) , |11)}. States resting on the red 
line connecting north and south pole have zero concurrence. 
States on the green equator are equivalent to Bell states thus 
maximally entangled. 



where s (t) is a random function. This is a classical 
noise model. To compute n (t) we need to average over 
a probability functional for s (t) , which we will specify 
more precisely below. Note that the manifold spanned 
by {|00),|11)} is decoupled from the manifold spanned 
by {|01), 1 10)} under the influence of this Hamiltonian. 
Thus if the initial density matrix lies in one of the two 
subspaces, the four-level problem decouples into two two- 
level problems and we can use the Bloch ball representa- 
tion to visualize the state space and entanglement evolu- 
tion. 

Take the initial state to be in span{|00), |11)} for ex- 
ample. The dynamical subspace D% is a 3-dimensional 
ball, as shown in Fig. [3] This makes it relatively easy to 
visualize the state and entanglement evolution. However, 
note that the center of D% is not the state n — 0. In fact, 
all the points on the z-axis belong to B$ because every 
neighborhood of any of these points contains points for 
which C > 0. 

The square roots of the non-zero eigenvalues of pp are 
Ai,2 = (Vl — r 2 cos 2 9 ± r sin6*)/2, where r and 6 are the 
spherical coordinates of the ball, and p is the spin-flipped 
density matrix, as in Eq. |H1 The concurrence is given by 



C = 



(15) 



H(t) 



1 



[,Ja A - <r B + s(t) go -Z + B o-*] 



(14) 



The maximally entangled states are on the equator and 
the separable states are on the z-axis, as seen in Fig. [3] 
The concurrence has azimuthal symmetry and is linear in 
the radial distance from the z-axis. The separable states 
in Z?3 form the 1-dimcnsional line D3 n S that connects 
the north and south poles of D3 . 

The key point is that D3 n S has a lower dimension 
than D3 itself. Now consider the possible trajectories 
with rioo on the z-axis. No function n(t) with continuous 
first derivative can have a finite time interval with C(t) = 
0. The trajectories either hit the z-axis at discrete time 
instants which puts them in category B, or approach the 
z-axis asymptotically which puts them in category A. 



Let us specify s (£) in more detail to demonstrate how 
those two qualitatively different behaviors are related to 
Markovianity. Qubit A sees a static field Bq while qubit 
B sees a fluctuating field s(t)g. All fields are in the in- 
direction. We will take the noise to be random telegraph 
noise (RTN): s(t) assumes value ±1 and switches between 
these two values at an average rate 7. RTN is widely 
observed in solid state systems |46l449| . 

For this dephasing noise model, the above-mentioned 
decoupling into two 2-dimensional subspaces occurs. In 
the 2 x 2 block labelled by {|00) , |11)} we find 



H(t) = -^-[a(t)g + B } 



wherer a z is the Pauli matrix in the subspace. 

This Hamiltonian can be solved exactly using a quasi- 
Hamiltonian method [32]. The time-dependent deco- 
herence problem can be mapped exactly to a time- 
independent problem where the two-value fluctuating 
field is described by a spin half particle. The quasi- 
Hamiltonian is given by 



H q = 



-ij + i-yn + L Z (B + r 3 g), 



where Ti are the Pauli matrices of the noise "particle" . 
L z is the SO(3) generator in the {|00) , |11)} space. 

The transfer matrix is given by 



T(t) = 



( T {t)cosBot £ T (t)smBot 6 

-( T (t)smB t ( T (t)cosB t 

1 



(16) 



where 



Cr(i) 



-7* 



cos (Tit) + j- sin (fit) 



Q= ^T^2 



(17) 
(18) 



is the dephasing function due to RTN and it describes 
the phase coherence in the x-y plane (l9j |. m(t) = 
since the dynamics is unital. Note (r(£) has qualitatively 
different behaviors in the g < 7 and g > 7 regions, as 
the trigonometric functions become hyperbolic functions 



Taking |$+) = (| 00) 
effective Bloch vector is 



|ll))/\/2 as initial state, the 



i(t) = 



(t (t) cos B t 

-CrW sin B t 





(19) 



The state trajectory is fully in the equatorial plane, as 
seen in Fig. @] The dephasing function £t(£) modulates 



the radial variation and the static field Bq provides pre- 
cession. 



m 



ICtWI, 



(20) 



and the dynamics is distance Markovian if g < 7 (£r (t) 
being monotonic). In this parameter region, Ct(£) can 
be approximated by exp(— g 2 /2 r y) and the dynamics is 
approximately Markovian as well. Thus we do not need 
to distinguish Markovian and distance Markovian in this 
model. In the Markovian case, the monotonicity of (t 
gives rise to spiral while in the non-Markovian case the 
state trajectory periodically spirals outwards with fre- 
quency £7. In both cases, the limiting state is the origin 
of the ball. 

The concurrence evolution is given by 

c(t) = Kt(*)|. 

Thus Markovian noise gives rise to entanglement evo- 
lutions in category A while non-Markovian noise gives 
rise to that in category B. Entanglement evolutions in 
the other two categories can not occur due to fact that 
dim (Dans') < dim(£> 3 ). 



C. Model D s 

In the previous section we saw that dynamical sub- 
spaces spanned by two computational basis states does 
not possess the property dimD n S = dimD. A natural 
question is whether simply increasing the number of basis 
states helps. This can be done by choosing a Hamilto- 
nian that connects only the triplet states in the original 
Hilbcrt space. One example would be 

H = J(t) a A -a B , 

as in Ref. [501, where the Heiscnberg coupling J it) has 
time dependence and is modeled as a classical random 
process. 

Note this Hamiltonian conserves the total angular mo- 
mentum of the two qubits. As a result, the triplet space 
spanned by {|00) , |11) , (|10) + |01)) /\/2} is decoupled 
from the singlet space. Thus the dynamical subspace 
D has its basis elements in su(3) if we choose the initial 
state to be in the triplet subspace. 

Using Gell-Mann matrices as the elements of su(3) al- 
gebra, the state space is a subset of a ball in R 8 (29| . 
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The rrii are linearly related to the previously defined rij 
The square roots of the eivenvalues of pp are 
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Markovian 
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FIG. 4. (Color online) (a) Entanglement evolution and (b), (c) state trajectories in the dynamical subspace Dz with RTN 
for Markovian and Non-Markovian evolutions. Concurrence is equal to the absolute value of Ct(4)- State trajectory in the 
equatorial plane of the effective Bloch representation where the only separable state is the origin 0. The initial state is the Bell 
state |$+). The concurrence is given by the radial distance. In both cases, g = 0.5, Bq = 1. 7 = 1 for the Markovian case and 
7 = 0.1 for the Non-Markovian case. 



Ai, 



mi + mi ± 



y/2 



2 — v3 ( v3m 3 + m 8 j + 3m 8 I v3: 



m 3 - mg 



^3,4 



(21) 



Thus the set D s n S of separable states is composed of 
two geometric objects: m§ = m^ = which is in R 6 and 



V3 ( v3rn3 + mg j + 3m 8 ( V3m 3 — mg ) =0 



which is in H. . In addition to a concurrence-zero 
hyperline, we have a C — hyperplanc in _D 8 . 
Hence dim (flg fl 5) < dim (D$) and this model only dis- 
plays entanglement evolutions in categories A and B. 

Introducing extra dimensions thus helps to form non- 
zero volume of D n S in D but a Hilbert space spanned 
by three of the four computational basis states is still 
not enough. D% and D% both avoid the region near the 
fully mixed state 0, where most separable states reside 
|5l| . This region is included in the dynamical subspaccs 
in the next two sections. 

Note that for D$ and D% , the symmetry of the Hamil- 
tonian and the specification of the initial conditions allow 
us to fully describe the dynamical subspace without ex- 
plicitly solving the dynamics. This feature can be seen 
in the more complicated models in the following sections 
as well: entanglement evolution categories, as a qualita- 
tive property of the system dynamics, can be determined 
from symmetry considerations of the model (dynamics 
plus initial condition), position of rioo in D PI S and the 
memory effect of the environment. 



D. Model YE 

Yu and Eberly considered a disentanglement process 
due to spontaneous emission for two two-level atoms in 
two cavities. In this case the decoherence clearly acts in- 
dependently on the two qubits. Nevertheless they found 
that ESD occurs for specific choices of initially entangled 
states [g. 

The decoherence process is formulated using the Kraus 
operators 

4 

p(t) = J^K„(t)p(0)Kl(t), (22) 

where K^t) satisfy J2u K l K » = I for a11 * 0112. For 
the atom-in-cavity model, they are explicitly given by 



where 



and 7 = exp(— Tt/2) and uj = y/l — exp(— Tt). 

It is possible to choose initial states such that the den- 
sity matrices have the following form for all t 

■a(t) 

b(t) z(t) 

z(t) c{t) 

d{t) 



K l =F 1 


®Fi 


K 2 = 


F X ®F 2 


(23) 


K 3 =F 2 <g> Fi 


Ki = 


F 2 ®F 2 , 


(24) 


Fi = 


7 0" 
1 


,F 2 = 


'0 0" 
u 


3 


(25) 



,«) = \ 



(26) 



with 



a(t) = K 2 a , (27) 

b(t) =c(t) = k + k(1-k)oo, (28) 

d(t) = l-o H-2(l-«) + (l -K) 2 a , (29) 

z{t) = k. (30) 

where k = exp(— Ft). Here the parameter ao determines 
the initial condition. 

The two-qubit entanglement is 



C(t) = - max {(),«/(*)}. 



(31) 



where /(£) = 1 - v/a [l - a + 2(1 - k) + (1 - n) 2 a }. 

In the polarization vector representation, the dynamics 
defined by Eq. [52] can be given explicitly by the transfer 
matrix T(t) and the translation vector rn(t) 



T(t) 



( « n 



w 



K 





K 





Ok 








K 


-K 


K 2 - 



' 





K K 2 / 



(32) 



and 



i{t) = [k- 1;0;0;k- 1;(k- l) 2 ]. 



(33) 



Here the coordinates are {niz, nxx, nyy, nzi,nzz}- 
We note that Eqs. [T2l and [TBI arc satisfied and this sponta- 
neous emission model is Markovian and distance Marko- 
vian. 

The non-zero components are 



niz(t) =n Z i(t) =-! + -(! + a Q )n (34) 



n X x(t) =nyy(*) = -K 

4 4 2 

^zz(i) =1 - r(l + Oo)« + o a o K • 



and the limiting state is 



■1,0,1), 



(35) 
(36) 

(37) 



where the coordinates are (njz, nxx , n zz)- 

This shows that the dynamical subspace Dy is a 3- 
dimensional section of M where the non-zero components 
are n IZ ,nzi,n X x,nYY, and n zz but also n IZ = n ZI 
and nxx = nyy, such that it can also be visualized in 
three dimensions. Interestingly, the limiting state rioo 
due to spontaneous emission is on the boundary of set S 
of separable states, and the purity of the state increases 
with time. 

Although we have so far fully solved the system dynam- 
ics, for the purpose of describing Dy , it is enough to know 
that {IZ, ZI, XX, YY, ZZ} are the basis of D Y and that 
the Kraus operators preserve the equalities niz = nzi, 
nxx = nyy from the initial conditions. Dy is then 




XX 



FIG. 5. (Color online) The dynamical subspace Dy is a tetra- 
hedron outlined by black lines. The set of separable states 
Dy n S form a hexahedron outlined by green lines. The ori- 
gin (fully mixed state n — 0) is denoted by a red filled dot. 
The limiting state ftoa in Ref. 0| is denoted by a blue circle. 



determined from the positivity condition of the density 
matrix. 

The positivity condition for the density matrix is given 
by 



0-2 > =► 2n IZ 



2n xx 



l zz 



<3 



a 3 > =*> 2nj z + 2n 
a 4 > =► AB > 



L xx 



i 2 zz < 1 + 2n zz {n 2 IZ - n 2 xx ) 



where A = (1 + n zz ) 2 - 4n 2 IZ and B = (1 - n Z zf - 
An xx . At t = 0, Ho = (2a Q - 1, 2, -l)/3. The positivity 
constraint a^ > gives rise to the range of the possible 
initial conditions, parametrized by ao G [0, 1]. 
The square roots of the eigenvalues of pp are 



A„ =X b = -VA 

A c =- |1 -n Z z + 2n X x\ 

Ad =-. |1 -n Z z - 2n X x\ 



(38) 
(39) 
(40) 



The ordering of the A's can change during the course of 
an evolution. When A a is the largest one finds C = 0, 
which is helpful in dtermining Bg. 

Dy is a tetrahedron with vertices at (0,-1,-1), 
(0,1,-1), (-1,0,1), and (1,0,1), as seen in Fig. M 
Dy n S is a hexahedron that shares some external areas 
with Dy. The 2-dimensional section of Dy with niz = 
is shown in Fig. [6l On the other hand, if the section is 
done with nxx = 0, we get a upside down triangle made 
of separable states, as shown in Fig. [7] 

Yu and Ebcrly showed that a sudden transition of the 
entanglement evolution from category A to category £ 
is possible as one tunes the physical parameter ao- This 





XX 



FIG. 8. (Color online) Skematic drawing of the different types 
of entanglement evolutions in the dynamical subspace Dy. 
The limiting state n M is denoted by a blue circle on the 
top right corner. Entanglement evolutions in category A are 
brown-colored while those in category £ axe black-colored. 



FIG. 6. (Color online) Cross section of the dynamical sub- 
space Dy with niz = and Dz with uxy — 0. 0,2 — is 
the ellipse. 03 = gives the parabola and the bottom of the 
isosceles triangle. 04 = sets the two sides of the triangle. 
The entangled region is shaded where the filled arrows denote 
increasing direction of the concurrence. The green dashed line 
is the boundary of entangled and separable states. The fully 
mixed state is denoted by a red dot. 




FIG. 7. (Color online) Cross section of the dynamical sub- 
space Dy with nxx = 0. 02 = is the ellipse. 03 = gives 
the parabola and the top of the isosceles triangle. 04 = sets 
the two sides of the triangle. The whole section is filled with 
separable states. 



phenomenon can be easily understood in our formalism, 
as seen in Fig. [3] The curvature of the state trajectories 
vary as the initial state changes. Thus there is a con- 
tinuous range of initial states parametrized by ao whose 
trajectories enter S fl Dy within finite amount of time 
and also a continuous range of initial states whose tra- 
jectories never enter Int(5). Note Fig. [5] is a schematic 



drawing since the true state trajectories are truly three 
dimensional. 

To be more quantitative, the transition between the 
category A and category £ behaviors in the YE model 
could be determined by examining the angle 600 between 
B$ and the tangent vector of the state trajectory in the 
long time limit. We denote the limiting tangent vector 
by ut(oo) and it is given by 



n T (oo) = (1 + a , l,4ao - 2). 



(41) 



The relevant B$ in the YE model is a plane passing 
defined by the following three points: n^, (0, 0, —1) and 
(0, 1/2,0). It is parametrized by 



(ft- noo) = 0, 



(42) 



where rh = (2, —2, l)/3 is the unit normal of the plane 
pointing into the separable region S H Dy- Note rh ■ 
71t(oo) = (6<2o — 2)/3 is proportional to cos #00 and its 
sign tells us whether the state trajectory approaches hoo 
from the separable region or the entangled region. Since 
< ao < 1, m • tit (00) falls in the range [—2/3,4/3] and 
both the ESD and HL behaviors are possible. 
The condition 



rh ■ tit(oo) = 0, 



(43) 



i.e., ao = 1/3, gives rise to the critical trajectory which 
approaches hoo along Bs- When rh ■ 7?t(oo) > 0, i.e., 
ao € [1/3,1], the state trajectory approaches hoc from 
the entangled region and we get entanglement evolutions 
in category A. These trajectories are represented by 
the brown curves in Fig. [5) On the other hand, when 
rh ■ hxioo) < 0, i.e., ao € [0,1/3), the state trajectory 
approaches hoc from the separable region and we get en- 
tanglement evolutions in category £. These trajectories 
are represented by the black curves in Fig. [5) 

The key point about the model YE is that the limiting 
state hoo G Bs'- it is on the boundary of the entangled 
and separable regions. That is why entanglement evolu- 
tions in both category A and category £ are possible. 
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E. Model ZJ 

Here we present a physically motivated dynamical sub- 
space where the dynamics satisfies the following condi- 
tions: 1) the two qubits are not interacting; 2) the noises 
on the qubits are not correlated; 3) the effect of dephas- 
ing and relaxation can be separated. This model shows 
all four categories of entanglement evolution. 

For this model the two-qubit dynamics can be decom- 
posed into single-qubit dynamics |19f . The extended two- 
qubit transfer matrix is T(t) = R ® R where 



R(t) 



n o 

(^(t) COS B()t t, \t) Bill JDQL 

-C(t)sinB < C(f)cosB i 
n 





C(t) sin B t 








-rit 



is the extended transfer matrix of individual qubits. The 
top left entry 1 describes the dynamics of / and is there 
only for notational convenience. Here £(£) describes de- 
phasing process, Ti is the longitudinal relaxation rate 
and Bq is static field in the z direction that causes Lar- 
mor precession. £(0) = 1 and C(°°) = if dephasing 
occurs. Note this dynamical description of decoherence 
is completely general as long as one can separate dephas- 
ing and relaxation channels. 

The dynamical subspace in this model is a specially 
parametrized 3-dimensional section of the full two-qubit 
state space M. Only the components nxx, tixy, nyx, 
nyy and nzz are non-zero and we further have con- 
straints nxx — tiyy and nxy = -iyx- We thus 
use nxx,XY,zz as independent parameters and the state 
space can be visualized in three dimensions. This dy- 
namical subspace has been previously considered in Ref. 
|42| and we will call it Dz- 

The fact that dim D < dim M relies on judicious choice 
of the initial states po ■ In Ref. [19| , more general initial 
states are considered such that Dz is expanded into a 
dynamical subspace with 7 elements in the su(4) algebra, 
i.e., {IX, XX, XY, YX, YY, ZI, ZZ}. 

Note n xx + n XY = R 2 is conserved in Dz- The posi- 
tivity of the density matrix requires 



a 2 >0 
a 3 >0 
a 4 >0 



2R 2 

z 

2R + n zz <l 



n% z < 3 
nzz < 1 — 2i? 2 , and n Z z > -1 



The concurrence is given by 



C = max \ 0, R 



1 + n Z z 



(44) 
(45) 
(46) 



(47) 



Separable states form a spindle shape on top and en- 
tangled states form a torus-like shape on bottom with 
triangular cross sections. A section along nxy = is 
shown in Fig. [6] 

We have thus fully described the entanglement topol- 
ogy of Dz and now we construct entanglement evolutions 



that induce Dz- A model similar to that of Eq. [Mlthat 
satisfies the three conditions is 



1 



H(t) = --[s(t)g-a B + B a*] 



(48) 



Note the RTN has both dephasing (g z ) and relaxation 
(g x ,gy) effects in this case. Situations with uxy = at 
intermediate qubit working point (with the presence of 
both dephasing and relaxation noise) have been consid- 
ered in Ref. 1421 • Here we choose the initial state to be 



the generalized Werner state [41 



u;*=r|$}<*| + 



1-r 



where 



|fc) = -(|00> + e*|ll» 

is a Bell state. 

The state trajectory is then given by 

n X x(t) =rcos(2B t)((t), 

n XY (t)=-rsm(2B t)((t), 

n Z z{t) =re- Tlt . 



(49) 



(50) 



(51) 
(52) 
(53) 



Similarly, if the Werner state derived from the Bell 

state 



|*> = i(|01>+e*|10» 

is used as initial state, the state trajectory is 



n X x{t) =rCr(t), 
nzz(t) = 



,-r-it 



(54) 



(55) 
(56) 



In both cases, the evolution of the concurrence is 

C(i) = max{0,r[|C(i)| -£(*)]} 

where £(£) = (1 — re~ Tlt )/2. Thus dephasing (relaxation) 
drives the state horizontally (vertically) towards the ori- 
gin. Note ESD occurs whenever £(i — > oo) ^ and the 
corresponding limiting state is n^ = |1{|. In the case 
of pure dephasing, T 1 = and the trajectory moves in 
the horizontal plane nzz = —T- The limiting state is 
«oc = (0, 0, — r). The only case where ESD does not hap- 
pen is the pure dephasing processes on the nzz = — 1 
plane since ftoo G Bs- This amounts to £(£) = and 
there is no abrupt cut-off. 

At the pure dephasing point, \n(i) — n^l = v2r|£(i)| 
and the dynamics is distance-Markovian if £(£) is mono- 
tonic. Since the typical Markovian dynamics can be char- 
acterized by £(£) with exponential decay forms, we do not 
distinguish Markovian and distance-Markovian for this 
model. Entanglement evolutions in categories £ and O 
are possible, depending on whether the RTN is Marko- 
vian or non-Markovian, as seen in Fig. [S] 
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FIG. 9. (Color online) (a) Entanglement evolution and (b), (c) state trajectories in Dz with pure dephasing noise. Only one 
qubit is subject to RTN and the other sees a static field Bo- State trajectories with initial state wf + is colored brown while 
those with initial state wf + are colored black. £(t) is given by Eq. 1171 and 7 = 0.5, g — 0.1, r = 0.5, Bo = 0.1. 



If the initial state is on the nzz = — 1 plane, by tun- 
ing the relative weight of dephasing and relaxation noise, 
the entanglement evolution can transit from category A 
to category £ if the noise is Markovian, or from cate- 
gory B to category O if the noise is non-Markovian (42| ■ 
These transitions can be interpreted as a transition of 
entanglement topology of the dynamical subspace as the 
physical parameter varies. At the pure dephasing point, 
Dz shrinks to the nzz = — 1 disc and ftoo € Bs- When 
this occurs we find A and B behavior. 



IV. EVOLUTION CATEGORIES 

Having given examples of models that produce the 
most important categories of entanglement evolution, we 
now provide a proof that the four categories essentially 
exhaust the possibilities. There are three cases. 

Case 1. ?ioo £ Int(S'), Image(n) and Bs are transversal 
everywhere. Here ft : [0, 00) — > M is considered as a 
map from the time domain to the state space. 

Consider a single trajectory n (t) . The 1-dimcnsional 
manifold Image (ft) C M and the 14-dimensional man- 
ifold Bs satisfy dim(M) = dim [Image (ft)] +dim(.Bs). 
Hence we can apply mod 2 intersection theory to the in- 
tersections of these two manifolds [53[ • First assume that 
Image(n) and (Bs) arc transversal everywhere. (We dis- 
cuss this assumption further below.) Since M is convex, 
we are guaranteed that there is a straight-line path from 
ft (0) to ftoa and that the actual ft (t) is homotopic to this 
path. The number of times that the straight-line path 
crosses Bs is one. The mod 2 intersection theorem then 
states that the number of times that ft (t) crosses Bs is 
odd (formally, this is referred to as the cardinality of the 
0-dimensional manifold ft" 1 (Bs))- This yields two cat- 
egories as we have defined them: card {ft" 1 (Bs)} = 1 
gives category £, while card {ft" 1 (Bs)} = 3,5,... gives 
category O. 

Case 2. ft^ e Bs, Image(n) and Bs are transversal 



everywhere. 

We are again assured of the existence of the straight- 
line path. However, this path can belong either to cate- 
gory A: the path does not intersect Bs and the approach 
to ftrx, is from M\S, or to category £ : the path intersects 
Bs once, and the approach to ft^ is from S. Because S 
is convex, multiple intersections of the straight-line path 
with Bs are not allowed. All other paths arc again homo- 
topic to the straight-line path. The intersection theorem 
implies that if a general path approaches ftx, from M\S, 
then it belongs to category A or O (in the variant where 
C (t) > for t > tf, for some finite tf), while if the path 
approaches ft^ from S, then it belongs to category A or 

(in the variant where C (t) = for t > tf). 

Case 3. Image (ft) and Bs are not transversal every- 
where. 

There is the possibility that Imagc(n) and Bs arc 
not transversal at discrete values of t. This gives rise 
to quadratic or higher-order zeros of C (i); the occur- 
rence of such points is not of great interest since they 
lie in a set of measure zero that is not related to phys- 
ical conditions. Of greater interest is the case when 
non-transversal points of intersection are due to restric- 
tions placed on the paths that stem from physical re- 
strictions placed on a model, e.g., a symmetry. For a 
given model we must then consider D, its dynamical 
subspace (a precise description of D is given in Ap- 
pendix EJ. We must replace M by D, S by D C\ S, 
and Bs by Bs H D. The dimensions of these spaces are 
model-specific. For mod 2 intersection theory to apply, 
we need dim (D) = dim [Image (ft)} + dim (Bs flO) = 

1 + [dim(S D D) - 1] = dim(5 (ID). If this holds, then 
the reasoning applied to cases 1 and 2 can be repeated. 
However, it is also possible to have dim(iS'n.D) < dim (D) 
- we have seen an example in Fig. [3] In this case the in- 
tersections of Image(n) and Bs<~)D are never transversal. 
A generic path ft (t) must either approach or "pierce" the 
low-dimensional space. This gives rise to A and B behav- 
iors, respectively. 
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dim(DnS) = dimD 


dim(D nS)< dim D 


n^o € Int(S) 


£ or O 


AorB 


iioo € Bs 


A, B,£orO 


AorB 



TABLE II. The relationships between entanglement evolution 
categories and details of the dynamics, such as the position 
of the limiting state iioo , etc. 



V. RESULTS AND CONCLUSIONS 

We now summarize the results exemplified in Sec. IIIII 
and formalized in Sec. IIVI by outlining a protocol to 
determine entanglement dynamics of any physical model. 

For any physical model, the possible entanglement evo- 
lution categories are determined by three factors: the 
limiting point n^, the dynamical subspacc D, and the 
presence or absence of distancc-Markovianity. Once these 
questions are resolved, the possible entanglement evolu- 
tions arc given by Table HU 

Tioo is the fixed point of the asymptotic dynamics. 
Since it is determined by the long-time limit of the dy- 
namics, it is usually fairly easy to compute. It may be 
a function of the parameters of the model. The cru- 
cial question about n^ is whether it belongs to Int(S'), 
the interior of S (Row 1 of Table HTf or Bs, the bound- 
ary of S (Row 2 of Table [TTJ) . D is the collection of all 
trajectories of the model. It depends on all parameters 
of the model, and on their range of variation. A pre- 
scription for computing D has been given in Sec. IIIII 
A. The crucial questions about D are its dimensional- 
ity dim(D) and the dimensionality of dim(D nS). If 
dim(D) = dim(D n S) then Column 1 of Table M ap- 
plies. If dim(D) < dim(DnS) then Column 2 of TableU 
applies. Finally, non-distance-Markovian evolution gives 
B and O, while distancc-Markovian evolution typically 
gives £ and A. 

Once the possible behaviors are known, one would like 
also to be able to compute the critical points or surfaces 
where transitions from one to another occur. 

The first type of transition occurs when ftoo does not 
move between Bs and Int(S'), i.e., we stay in one row 
of the table. One cannot move between columns within 
a given model since D, by definition, represents all pos- 
sible states of a model. Thus we stay in one box of the 
table. £ < — > O and A < — > B are the possible transitions 
and the critical points represent a change from distancc- 
Markovianity to non-distance-Markovianity. Such tran- 
sitions may be subtle and difficult to calculate, and an 
explicit equation for the critical surface cannot be given. 
A quite simple example where the calculation can be done 
analytically has been given in Sec. IIIII B. 

The second type of transition happens when ft oo moves 
between B s and Int(S) and dim(D) = dim(DnS) (Col- 
umn 1). There will be a specific point at which this oc- 
curs. This point is determined by examining the asymp- 
totic dynamics of the model as a function of the parame- 
ters. If the parameters are collectively given by a set {gi}, 



then the equation for the critical surface is ftoo(gi) £ Bs- 
The transitions are £ < — > A and O < — > B for distance- 
Markovian and non-distance-Markovian evolution. An 
example has been given in Sec. IIIII E and Ref. J42J . 

The third type of transition occurs in Column 1, Row 
2 of the table. In this case n^ does not change, but the 
direction of approach to n^ changes as a function of the 
parameters. In this case the equation of the critical sur- 
face is ?T/r(oo) G tang(i?s), where nr(oo) is the limiting 
tangent vector of the state trajectory and tang(i?s) is 
the tangent space of the boundary of S. The transitions 
are £ < — > A and O < — > B for distancc-Markovian and 
non-distance-Markovian evolution. An example has been 
given in Sec. IIIII D. 

"Tricritical" points can also occur when more than one 
critical condition is satisfied. 



Appendix A: two dimensional sections 

With our conventions for the two-qubit su(4) algebra, 
the structure constants can be calculated from 



Jijk 
Qijk — 



1 



2% x 4 
1 



Tr \jM,tlj]Hk 



2x4 



Tr [fM,Hj] + fik, 



(Al) 
(A2) 



Here we denote [, ]+ as the anti-commutation operation. 
Most of them are zero and the non-zero structure con- 
stants are 

1 = ^1,4,5 = ^1,8,9 = ^1,12,13 = ^2,4,6 = ^2,8,10 

=^2,12,14 = G?3,4,7 = ^3,8,11 = ^3,12,15 = — ^5, 10,15 
— ^5,11,14 = <^6,9,15 = -^6,11, 13 = — ^7,9,14 = "7,10,13 
=/l,2,3 = /l,6,7 = /l,10,ll = /l,14,15 = — /2,5,7 
= — /2,9,11 = — /2, 13, 1.5 = J3,5,6 = J3,9,10 = /3, 13,14 
=/4,8,12 = /4,9,13 = /4,10,14 = /4,11,15 = /5,8,13 
=/5,9,12 = /6,8,14 = fe, 10,12 = f-,8,15 = h, 11,12- 

Suppose the two non-zero components of a two- 
dimensional section assume numerical values x and y. 
It can be shown that 



Trp 2 =l(l + i? 



Tr p 6 



■so 

1 



3i? 2 ) 



(A3) 

(A4) 



Tr p =— [1 + 6-R +R + ninjnen m dijkdki m }, (A5) 

where R 2 — x 2 + y 2 . Note the generators p,i as defined 
in Eq. [3] are elements of the Pauli group thus they either 
commute or anti-commute. This is the origin of the two 
geometrically different 2-dimcnsional sections in M. 

If the two Hi's anti-commute, the dijk related term 
vanishes and the positivity constraints is solely on 03, 



a 3 >0 



ir < 1. 



(A6) 



13 



This is the disc section. On the other hand, if the two 
/ij's commute, the d^ related term gives 4x 2 y 2 . The 
positivity constraint is on 04 

a 4 > => R 4 - 2R 2 - 2.tV + 1 > 0. (A7) 

This gives rise to the square section. 

Appendix B: volume estimation 

The sets M and S are spheroids of varying radius and 
the set of entangled states M\S is a spheroidal shell of 
varying thickness. To get some idea of the volumes of M, 
S and M\S and their radial profiles we have performed 
Monte Carlo simulations on the sphere of radius a/3 in 
K 15 . The results are shown in Fig. [TUJ where we give 
the probabilities that a state n is physical (positive den- 
sity matrix) and that it is separable, as a function of \n\, 
leading to a radial probability distribution function for 
the physical states and the separable states. The distri- 
butions fall off fairly abruptly, which allows us to identify 
a rough "effective" radius of about a/3/ 2 for the physi- 
cal states and about 0.45 x a/3 for the separable states. 
Plotted in this way, the entangled states form a rather 
thin shell. This is not the only way to display the data, 
of course. If the volume densities were to be plotted in- 
stead, the function shown would be multiplied by \n\ , 
and the density of separable states would show a sharp 
peak. 

To estimate the volume of separable and physical states 
inside the 15 dimensional ball with radius a/3, we first 
perform simple sampling Monte Carlo to get the ra- 
dial distribution function p(r) at various r values. Here 
r = \n\. 10 11 sample points are generated at each 
Tk = \/3 x 0.02A:. The accuracy of each p(rk) is beyond 
five digits. Uniformly distributed random numbers on a 
14 dimensional sphere is generated using the algorithms 
inRef. II. 

Take p(r) to be the radial distribution for the physical 
states for example. To get Vph ys from p(r), we perform 
numerical integration using Simpson's rule. 



1/ r l W Jo 3 p{r)r u dr 
1W \p(r)\ = - y — 7= 



V B 



=15Vb 



p(f)f dr. 



(Bl) 
(B2) 



where f = r/v3 is the normalized radius, i.e., |n|/i/3. 
Vb it the volume of a 15-dimensional ball with radius a/3 



V b = V(r = a/3, a = 15) 

W 2 



r( s /2 + i) 



1444.905 



(B3) 
(B4) 
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FIG. 10. Radial probability distribution functions of physical 
states (blue circle) and separable states (green square) in the 
ball of radius y3 in R 15 . For any fixed radius |n|, the vertical 
span of the white, yellow and cyan regions denotes the prob- 
ability of a randomly generated n being separable, entangled 
and unphysical. The volume of S and M can be calculated 
from Eq. lB2l 



From the radial probability distributions of Fig. [TUl 
we get the volume of separable states and physical states 
in the 15-dimcnsional polarization vector representation 



t'scp 


=0.008971 ± 0.000005 


(B5) 


Vphys 


=0.03700 ±0.00001 


(B6) 


^scp 


=0.2424 ±0.0002. 


(B7) 



v, 



phys 



The separable states only account for about a quarter 
of the state space in our volume measure, which is the 
Euclidean measure on M. This measure has the virtue 
that the distance between two different values of n cor- 
responds to a physical distance between 15 "pointers" if 
a complete measurement of all 15 observables is made. 
The volume in this measure is different from that of Rcf . 
[39j , where eigenvalues of the density matrix were used 
to discriminate states. With that approach the separable 
states occupy more than half the volume of the full state 
space. 

A probability distribution function P (C) for various 
\n\ can be computed in a similar way. This is shown 
in Fig. 1111 Again, appreciable entanglement is concen- 
trated at a fairly specific radius. At \ft\ = a/3, i.e., for all 
pure states, the measure of separable states is zero. This 
is because pure separable states are parametrized by 4 
parameters while pure states require 6 parameters. As 
stated above, S contains a ball of radius \n\ = 1/a/3 ~ 
0.57735. The numerical study indicates that the maxi- 
mum radius is at about \n\ ~ 0.5782. Although there is 
some angular variation in the function C(n), the magni- 
tude of n is a good rough indicator for the entanglement. 
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P(C) 



P(C) 




Appendix C: dynamical subspace D 

In the text we denned the dynamical subpace D infor- 
mally as the set of all trajectories allowed by the set of 
initial conditions and the equations of motion. Here we 
give a more formal definition. 

The authors of a model generally specify n parameters 
in the total Hamiltonian for the coupled system and bath. 
These parameters vary over some ranges, gi : mi n < g\ < 
Si, max, etc., giving a set G = {#;}■ Furthermore, a set 
of initial conditions is also given by a similar set H . For 
example, in the YE model, G — {T} with < T < oo, 
and H = {ao} with < ao < 1- In the literature, G 
and H are almost always manifolds (or manifolds with 
boundaries), as they are in this example. Finally T = {t} 
with < t < oo is always a 1-manifold. The equations of 
motion then define a smooth map N : Gx H xT — > M, 
where M is the state space. The dynamical subspace D 
is then defined asD = Image(-/V). 
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